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Abstract 

We present the first study to examine the validity of the relativistic impulse approximation 
(RIA) for describing elastic proton-nucleus scattering at incident laboratory kinetic energies lower 
than 200 MeV. For simplicity we choose a 208 Pb target, which is a spin-saturated spherical nucleus 
for which reliable nuclear structure models exist. Microscopic scalar and vector optical potentials 
are generated by folding invariant scalar and vector scattering nucleon-nucleon (NN) amplitudes, 
based on our recently developed relativistic meson-exchange model, with Lorentz scalar and vector 
densities resulting from the accurately calibrated PK1 relativistic mean field model of nuclear 
structure. It is seen that phenomenological Pauli blocking (PB) effects and density-dependent 
corrections to crN and wN meson-nucleon coupling constants modify the RIA microscopic scalar 
and vector optical potentials so as to provide a consistent and quantitative description of all elastic 
scattering observables, namely total reaction cross sections, differential cross sections, analyzing 
powers and spin rotation functions. In particular, the effect of PB becomes more significant at 
energies lower than 200 MeV, whereas phenomenological density-dependent corrections to the NN 
interaction also play an increasingly important role at energies lower than 100 MeV. 

PACS numbers: 24.10.-i,24.10.Ht,24.10.Jv,24.70.+s, 25.40. Cm,25.60.Bx 
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I. INTRODUCTION 



The relativistic impulse approximation (RIA) provides an excellent quantitative descrip- 
tion of complete sets of elastic proton scattering observables from various spin-saturated 
spherical nuclei at incident energies ranging from 200 to 400 MeV The latter represents 
the theoretical framework for generating complex microscopic optical potentials for solving 
the elastic scattering Dirac equation. The reliability of the RIA has also been demonstrated 
by the fact that, for the above energies, predictions of elastic proton scattering observables 
are very similar to the corresponding results based on the highly successful global Dirac 
phenomenological optical potentials, which have been calibrated to provide excellent quan- 
titative predictions of elastic proton scattering observables from stable nuclei ranging from 
12 C to 208 Pb and for incident energies between 20 and 1040 MeV 3|. 

The Melbourne group have already developed a highly predictive microscopic Schrodinger 
model for describing elastic scattering at energies between 20 and 800 MeV jjj, 4\. Here, 
however, we focus on the RIA which is based on the relativistic Dirac equation, and is 
more attractive in the sense that the microscopic scalar and vector optical potentials - 
and consequently the corresponding Schrodinger-equivalent central and spin-orbit optical 
potentials 1|- are directly related to the Lorentz properties of the mesons mediating the 
strong nuclear force, a connection lacking in nonrelativistic models. 

Recently we developed an energy- dependent Lorentz covariant parameterization of the 
on-shell NN scattering matrix at incident laboratory kinetic energies ranging from 40 to 
300 MeV 5j. In this paper, we employ the latter meson-exchange model to generate rel- 
ativistic microscopic scalar and vector folding optical potentials in order to systematically 
study the predictive power of the RIA for describing elastic proton-nucleus scattering observ- 
ables at energies lower than 200 MeV. At these low energies of interest, multiple scattering 
effects [6j, medium modifications of the NN interaction 3j, and Pauli blocking 

HQ 

con- 
tributions become increasingly important, and hence the validity of the RIA needs to be 
investigated at these low energies. 

In Sec. (jni), we briefly describe the RIA as well as the relevant input pertaining to our 
application thereof. Results are presented in Sec. (lllip . and we summarize and conclude in 

sec. am. 
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II. FORMALISM 



A. Relativistic Impulse Approximation (RIA) 



A comprehensive description of the RIA, as well as the corresponding computer codes 
and associated numerical details can be found in Refs. [l| and jg|. Here, we just give a brief 
introduction of it. The Dirac optical potential is given by 

A ■ A 
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where p is the magnitude of the three-momentum of the projectile in the nucleon-nucleus 
center-of-mass frame, |\&) is the A-particle ground state, and T is the nucleon-nucleon 
scattering operator. In the original RIA, T is chosen as: 



T(q,E) = Y,F L (q,E)Xf iy X 



L 

(2) 
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where scattering amplitudes F L (q,E) are complex functions of the momentum transfer q 
and laboratory energy E, which can be obtained by the relativistic Horowitz-Love-Franey 
model (see section Hi Bp . and Xh. stand for five Dirac matrices for the projectile or target 
nucleon. 

Assuming the ground state is a Hartree product of single-particle 4-component wave 
functions U a (r), the action of the optical potential on the incident wave \U ), projected onto 
coordinate space, can be written as: 

(v\U opt (q,E)\U (r)) = -^f £ [ f d\'p L (r')t L D (W - r|; E)] X L U (r) 

£ J d 3 r'p L (v\v)t L x (\v'-v\;E) 
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where 



d 3 q 
(2^ 



t L D (q,E)e-^ 



(4) 



with tp(q,E) = (iM 2 /2E c k c )F^(q) and similarly for the exchange pieces t x (Q,E). The 
nuclear densities are defined as: 



p L (r', r) = J2 ^(r')A^(r), /(r) = p L (r, r) 



(5) 
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Here the prime on the occupied states means that one sums over target protons when the 
density is to be used with pp amplitudes and over target neutrons when it is to be used with 
pn amplitudes. 

The first term in Eq. contains a multiplicative factor that defines the direct optical 
potential: 

UE(r,E) = -^f J d 3 r'p L (r')t L D (\r'-v\;E) (6) 
Adopting a local-density approximation, the second term can be localized to give: 

Uk(r,E) = -^ J dV/(r',r)t£(|r' - v\; E)j Q (p\v' - r|) (7) 

where jo is a spherical Bessel function. For spin-zero nucleus, the non-zero densities are only 
baryon, scalar and a tensor term associated with er°\ Hence, the RIA optical potential can 
be written as 

U opt (r; E) = U s {r- E) + j°U v (r; E) - 2ia ■ rU T (r; E) (8) 

where 

U L {r-E) = U L D {r-E) + U^r-E) (9) 

These optical potentials serve as input to solve the Dirac equation for elastic proton scat- 
tering so as to generate the relevant partial-wave scattering phase shifts for computing the 
scattering observables, namely the total reaction cross section or, the differential cross sec- 
tion da/dQ, analyzing power A y and spin rotation function Q. Similar to Ref. [lj], we resolve 
ambiguities in the form of the relativistic NN scattering amplitudes by employing pseudovec- 
tor coupling for the 7rNN vertex. In Ref. l|, Horowitz and Murdock have shown that the 
tensor potential has a negligible effect on all the observables and nuclei of interest. In a 
similar vain and for comparison to the latter, we also neglect the small tensor term. 

B. Relativistic Horowitz-Love-Praney model (HLF) 

The original Love-Franey model and its relativistic version, i.e., Horowitz-Love- Franey 



10 



111 ] and references therein. Here, 



model, have been described in detail in Ref. 
we briefly allude to list some important formula and the fitting procedure at lower energy. 
Essentially HLF model parameterizes the complex relativistic amplitudes F L (q, E) in terms 
of a set of = 10 meson exchanges in first-order Born approximation, such that both direct 
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and exchange NN (tree-level) diagrams are considered separately, that is: 

iM 2 

F L (q,E) = ^[F£(q)+F^Q)}, (10) 



where 

N 



FE(Q) = J2 S ^)(ri-r 2 ) T ^f(q) (11) 

i=l 

N 

F k(Q) = (-i) TNN J2 c mAn-T2) Ti nQ) • (12) 



i=i 



Here T« = (0,1) denotes the isospin of the i meson, Tnn refers to the total isospin of the 
wo-nucleon system, (fi ■ r 2 ) Ti is 1 or -3 for different Tnn and T i; Clu),l is the Fierz matrix 



m, and 



f ( x ) = "T^l 1 + TIT 2 - + S)" 2 ( 13 ) 

x 2 + mf Af x 2 + m 2X A 2 ' v 7 

where x represents the magnitude of either the direct three-momentum transfer q or the 
exchange-momentum transfer Q, (gf,gf), (rrii,rhi), and (Aj,Aj) are the real and imaginary 
parts of the coupling constant, mass, and cutoff parameter for the i meson. These param- 
eters are obtained by fitting the theoretical amplitudes with the values extracted from the 
NN scattering data. 

Recently we have developed a set of parameter at incident laboratory kinetic energies 
ranging from 40 to 300 MeV Q|, in which the coupling constant is set as energy dependent, 
namely 

g\E) = gl\l + a g (e^ - 1)] (14) 

where 

T rcl = (15) 

with To = 200 MeV, T rc i is positive in the 50 to 200 MeV energy range of interest, and g^, 
a g and ar are dimensionless parameters extracted by fitting to the relevant data. 



C. Relativistic Mean Field theory (RMF) 



For extracting the relevant scalar and vector proton and neutron densities to be folded 
with the HLF NN scattering matric, we employ the so-called PK1 Lagrangian density 13], 
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namely 



C 



vfd^ — M — g^cr — g^o^ - ^7 M f ■ p M - e7^ 

1 r\ r\n 1 2 2 1 3 1 4 

-d^ad^a - -m a a - -g 2 a - 7 g 3 a 



2 

4 ^ 



(16) 



where the parameters are listed in Table [H and the field tensors for the vector mesons and 
the photon are respectively defined as, 
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TABLE I: The PK1 effective interaction. The masses (in MeV), meson-nucleon couplings, and 
nonlinear coefficients are listed. 
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M p m a 
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92 


93 


C3 


939.5731 


938.2796 514.0891 


784.254 763 10.3222 


13.0131 


4.5297 


-8.1688 


-9.9976 


55.636 



Under the variation of the Langrangian density with respect to the different fields, one 
obtains the equations of motion for the nucleon and different mesons via a self-consistent 
procedure. To obtain the parameters, the masses of spherical nuclei 16 0, 40 Ca, 48 Ca, 56 Ni, 
68 Ni, 90 Zr, 116 Sn, 132 Sn, 194 Pb, and 208 Pb, the compression modulus K, the baryonic density 
at saturation p sat , and the asymmetry energy J of nuclear matter are fitted to give the min- 
imum by the Levenberg-Marquardt method. The relative deviations are 0.0102 and 0.0178 
for the masses of 19 spherical nuclei and charge radii of 13 spherical nuclei, respectively, 
which are similar or better than the other commonly used RMF effective interactions, such 
as the NL1 jy), PL-40 15], NL-SH 16], TM1 Q, NL3 Q, TW99 Q, and DD-ME1 Q. 
Compared to these effective interactions, the PK1 parametrization represents an improve- 
ment in the sense that center-of-mass corrections are included microscopically and it also 
provides a unified description of the properties of stable and unstable nuclei over a broad 
mass (16 < A < 214) range. 
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At last, to minimize uncertainties associated with nuclear structure input, we focus on 
proton scattering from 208 Pb which is a spin-zero spherical nucleus for which the PK1 ef- 
fective interaction can give excellent description. In Fig. [T], the scalar and vector density 
distributions for proton and neutron of 208 Pb are shown. Another advantage for considering 
such a heavy nucleus as 208 Pb, is that at the low energies of interest, recoil corrections to 
the Dirac scattering equation are expected to be small and hence can be neglected. 
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FIG. 1: (Color online) Scalar, and vector (baryon) density (in units of fm -3 ) distribution in radial 
direction r (in fm) for proton and neutron of 208 Pb, respectively. The proton and neutron densities 
are denoted by solid and dashed curves, respectively. The dash dotted curve corresponds to the 
proton scalar density, and dotted curve indicates the neutron scalar density. The PK1 effective 
interaction is used. 
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III. RESULTS 



Guided by the availability of experimental data at the lower, middle and higher-energy 
regions of the 50 - 200 MeV range, we focus on incident energies of 65 21|, 98 22j, 121 



231 ] and 200 MeV 24j, |25j, and examine the validity of our RIA predictions for describing 
da/dQ, A y and Q from a 208 Pb target nucleus as the energy is systematically lowered. Note 
that at approximately 100 MeV, there exists no complete set of experimental data for one 
specific energy, and hence we compare our RIA calculations (in the mid-energy range) to 
data for da/dQ and spin observables, A y and Q, at 121 MeV and 98 MeV, respectively. 

In Figs. [21 [31 and HI the RIA calculations are denoted by black dashed curves (color 
red online) and for reference we also display results based on the energy-dependent mass- 
independent (EDAI) global Dirac optical potentials (GOP) j^J indicated by the black dotted 
curves (color olive online). We consider center-of-mass scattering angles ranging up to a 
maximum value corresponding to a three-momentum transfer of about q ~ 2.5 fm -1 for 
which first-order microscopic nonrelativistic models and relativistic global optical potentials 



have been shown to be valid 



26]. As expected, the global Dirac optical potentials provide 



an excellent description of all the observables for the entire energy range of interest. At 
200 MeV, the RIA provides a satisfactory description of da/dQ and A y for angles ranging 
up to a value corresponding to the third minimum in da/dQ. Although the RIA provides a 
good description of Q from 20 to 40 degrees, it fails to describe the small angle behavior. At 
lower energies, the RIA provides a satisfactory description of da/dQ up to approximately 
the second minimum, but fails to describe A y and Q. 

For elastic proton scattering from a heavy nucleus such as 208 Pb, Murdock and Horowitz 
[l| have demonstrated that Pauli Blocking (PB) corrections to the RIA optical potentials 
play an increasingly significant role at lower energies. Essentially PB represents the mecha- 
nism that prevents nucleons in the nuclear medium from scattering to occupied intermediate 
states 2?|]. We now proceed to study the effect of PB corrections for our low energy range of 
interest. Following the procedure outlined in Ref. [lj], we include PB corrections to the RIA 
microscopic optical potentials via the so-called energy-dependent PB coefficients as,y(Ti a b), 
that is, 



1 - a s ,v(Ti 



lab J 



pB{r) 
Po 



2/3' 



U s ' v (r,T Hb ) 



(18) 
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where Ti a b denotes the incident laboratory kinetic energy, ps(r) is the local RMF baryon 
density of the target nucleus, po = 0.1934 fm~ 3 , and U s ' v (r, Xi a b) represents the uncorrected 
RIA scalar S and vector V optical potentials. Within the context of a relativistic Dirac- 
Brueckner approach, values for the PB correction factors a,s,v{Tiab) have been extracted at 
discrete energies of 135, 200, 300 and 400 MeV [jj and also at discrete momenta ranging 
from about 0.46 to 0.75 GeV/c [28]. Consistent with the results reported in Ref. [l], in Fig. [2] 
we confirm that PB corrections (black solid curve) provide an improved description of the 
200 MeV data, the effect been most pronounced for A y and Q at small angles. Due to the 
enhanced sensitivity exhibited by spin observables, we now focus on the influence of PB on A y 
at a lower energy of 98 MeV, where PB corrections are expected to be more significant: note 
that there are no experimental data for Q at 98 MeV. The lack of published PB coefficients 
at this energy necessitates linear extrapolation of the published Dirac-Brueckner values in 
Ref. [lj], which is shown in Fig. [5j RIA predictions with and without PB corrections are 
denoted by the gray solid (color cyan online) and black dashed (color red online) curves 
in Fig. [6] (a), respectively. One clearly observes that the extrapolated Dirac-Brueckner PB 
corrections cannot account for the A y data. 

Due to the latter failure, we now proceed to study PB via a phenomenological approach. 
In particular, we study to which extent phenomenologically extracted values of asy{T\ a \^) 
can provide a consistent description of complete sets of scattering observables at the specific 
energies of interest. Starting with the lower momentum values of as,v (Tiab) published in 
table 2 [po, P = 0.46 GeV/c] of Ref. 28[, we systematically vary the real and imaginary 
parts to give one GOP that provides the best fit-by-eye of da/dQ, A y and Q for elastic 
proton scattering from 208 Pb at five discrete energies of 65, 100, 135, 170 and 200 MeV, 
spanning the range of interest. In fact, this is enough for the purpose of this paper: to check 
the validity of the RIA at lower energies. Our phenomenological PB factors are listed in 
Table HT1 together with the corresponding published Dirac-Brueckner values from Ref. [l]. 

We now consider the effect of these phenomenological corrections on A y at 98 MeV, 
where the extrapolated Dirac-Brueckner corrections were shown to fail. In Fig. [6] (a) one 
clearly sees that our phenomenological PB corrections (black solid curve) provide a more 
satisfactory description of A y compared to the Dirac-Brueckner PB corrections (gray solid 
curve). Focusing on the optical potentials, in Fig. [6] (b) we observe that the main effect of 
both Dirac-Brueckner (grey solid curves) and phenomenological (black solid curves) PB is 
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TABLE II: Phenomenological Pauli blocking correction factors a(Ti a b) for real and imaginary RIA 
microscopic scalar S and vector V optical potentials as various energies T^. Where available, the 
corresponding published Dirac-Brueckner values [l] are indicated in square brackets. 



Scalar Vector 

Energy [MeV] Real Imaginary Real Imaginary 

65 0.020 0.68 0.14 0.85 

100 0.007 0.60 0.11 0.73 

135 0.005 [0.00377] 0.53 [0.108825] 0.09 [0.08403] 0.60 [0.24535] 

170 0.008 0.42 0.07 0.50 

200 0.010 [-0.0078] 0.35 [0.098] 0.0605 [0.0605] 0.42 [0.207] 

to decrease the strengths of the real and imaginary parts of the scalar and vector optical 
potentials, the effect being most pronounced for the imaginary parts. Furthermore, our 
phenomenological corrections have a larger effect than the corresponding Dirac-Brueckner 
corrections. Also note that the PB-corrected optical potentials differ in strength from the 
corresponding global optical potentials (black dotted curve), which suggests that the delicate 
interplay between the relative strengths and signs of the scalar and vector optical potentials, 
rather than absolute magnitudes, is responsible for providing an accurate description of the 
diffractive behavior of all scattering observables. 

In Figs. El 0, and HI we now display our theoretical predictions including the above- 
mentioned phenomenological PB corrections (black solid curves) over the entire range of 
interest. Although not shown, both Dirac-Brueckner and phenomenological PB corrections 
provide the same quantitative improvement for all observables at 200 MeV. At 121 MeV, PB 
corrections provide an excellent description of the first few maxima and minima for da/dQ. 
There are no data for Q at 98 MeV, but we note that PB corrections give results very similar 
to the corresponding global Dirac optical potentials. At the lower energy of 65 MeV, PB 
corrections invoke a larger effect and provide a satisfactory description of da/dfl and A y up 
to about 50 degrees, but predictions deteriorate at larger angles. Note, however, that at 
65 MeV, PB corrections fail to quantitatively describe the first two extrema for A y , as well 
as fail to reproduce Q over the entire angular range. The 65 MeV results emphasize the 
important fact that optical potentials can be deemed reliable only if complete sets of elastic 
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scattering observables, as opposed to only one or two observables, are accurately described. 

To further validate our phenomenological procedure for including PB corrections, in Fig. [7] 
we now consider a different spin-zero spherical target nucleus, namely 40 Ca, at an incident 
proton energy of 152 MeV: the meaning of the various curves is the same as those in Fig. 
Compared to the observables without PB (black dashed curves) , and consistent with our con- 
clusions for 208 Pb, the corrected values (black solid curves) provide an improved description 
of all observables and also agree with the GOP results. 

Although PB systematically improves the RIA predictions at lower energies, these cor- 
rections do not provide a quantitative description of complete sets of observables, and hence 
it is clear that other effects also begin to play an increasingly significant role as the energy 
is lowered. The question now arises as to whether one can incorporate additional correc- 



tions so as to consistent 
RIA. Various authors 



y 



30 



improve these low energy results within the framework of the 



31| have stressed that importance of nuclear medium modifica- 



tions to the NN interaction at lower energies. In particular, a number of theoretical models 
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32 



33 



34j predict density-dependent corrections to meson-nucleon coupling constants as 



well as nucleon- and meson-masses in normal nuclear matter. We now adopt a phenomeno- 
logical approach to investigate to what extent the renormalization of certain meson-nucleon 
coupling constants can provide a systematic improvement of the data at 65 MeV. Since 
scalar and vector RIA optical potentials are dominated by contributions from a- and uj- 
meson exchange to the NN scattering amplitudes, we vary the values of the real crN and 
u;N coupling constants, g 2 a and g 2 , so as to provide the best consistent description of elas- 
tic proton-nucleus scattering. The result of this procedure is that the value of g\ changes 
from -7.5701 to -11.1734 and g 2 changes from 7.4511 to 12.9704: these new values of the 
couplings correspond to a g = 0.0271 and cit = 3.8364 for the a meson and to a g = 0.1022 
and Or = 1.9527 for the uj meson [see Eq. (I14p in this paper, also Eq. (21) and Table II in 
Ref. [5]]. The effect of renormalizing these couplings is illustrated by comparing the original 
results (black solid curves) to the renormalized values (gray solid curve) in Fig. HI The most 
pronounced improvement is observed for the minima of da / dQ at large scattering angles and 
spin observables at small scattering angles. At the level of the optical potentials, the main 
effect of the modified coupling constants is to increase the strengths of the real parts of the 
scalar and vector potentials at lower energies such that the energy- dependent trend these 
strengths is qualitatively similar to the corresponding GOP results, as illustrated in Fig. [HJ 
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Furthermore, these improved predictions arising from renormalized coupling constants are 
consistent with other nuclear reaction- and structure-studies focusing on density-dependent 
corrections to the NN interaction 



13 



30 



31 



351 ] . Note, however, that the best descrip- 



tion of A y and Q up to angles of about 40 degrees, is obtained by including both PB and 
density-dependent corrections (gray solid curve). 

To further test the consistency of our scheme for including both PB and density-dependent 
corrections, we now focus on elastic proton scattering at an energy as low as 30 MeV where 
these effects are expected to be important. In Fig. [9], it is gratifying to observe that both 
phenomenological corrections (gray solid curves) contribute significantly toward providing 
an excellent description of the da/dQ and A y data 36|, |37J] over the entire angular region. 
Although there are no data for Q at this energy, we note that our predictions are similar to 
the GOP results which have constrained to fit elastic scattering observables at energies as 
low as 20 MeV Q]. The above mentioned results inspire confidence in the predictive power 
of our phenomenological corrections to the RIA. 

Finally, we check the validity of the RIA model for describing total reaction cross section 
o~r data for proton scattering from 208 Pb for energies lower than 200 MeV. This integral 
quantity serves as crucial input in designing accelerator-driven-systems for transmutation 
of radioactive waste. However, due to the shortage of relevant <tr data, simulations are 



forced to employ model-dependent values as input [38J. Furthermore, cxr values also form 
an essential ingredient in studies of nuclearsynthesis and stellar evolution. Consequently, 
the need for developing predictive models of elastic proton scattering should be evident. In 
Fig. [Til] we observe that compared to the uncorrected RIA predictions, denoted by the black 
dashed line (color red online), the inclusion of both PB and density-dependent corrections, 



indicated by the gray solid line (color cyan onl ine), pr ovides a satis: 



39 



40 



41 



42 



43 



44 



actory description of 



45 



461 ] . Our results 



261 ] as well as the 



the data: the experimental data are from Refs. 

are consistent with corresponding nonrelativistic microscopic predictions 
GOP values (black dotted curve), over the entire energy range. In particular, both RIA and 
nonrelativistic models fail to describe the 70 MeV data: a point of concern, however, is the 



unsystematic trend exhibited by the data in this region, as already mentioned in Ref. 



26]. 
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IV. SUMMARY AND CONCLUSIONS 



We present the first study to examine the validity of the relativistic impulse approximation 
(RIA) for describing elastic proton-nucleus scattering at incident laboratory kinetic energies 
lower than 200 MeV. For simplicity we considered a 208 Pb target, which is a spin-zero 
spherical nucleus for which accurately calibrated relativistic mean field nuclear structure 
models exist, and for which recoil corrections to the Dirac scattering equation are expected 
to be negligible. Microscopic scalar and vector optical potentials are generated by folding 
our recently developed relativistic meson-exchange model of the NN scattering matrix with 
appropriate Lorentz densities arising from the PK1 Lagrangian density. We have established 
that phenomenological Pauli blocking effects and density-dependent corrections to the erN 
and u>N meson-nucleon coupling constants modify the RIA microscopic scalar and vector 
optical potentials so as to provide a consistent and quantitative description of all elastic 
scattering observables (<jr, da/dfl, A y and Q) at energies ranging from 30 to 200 MeV. 
In particular, the effect of PB becomes more significant at energies lower than 200 MeV, 
whereas phenomenological density- dependent corrections to the NN interaction also play an 
increasingly important role at energies lower than 100 MeV. Note that although our initial 
study has been of a phenomenological nature, our results clearly indicate the importance of 
including PB and density-dependent corrections to the microscopic RIA model for providing 
consistent and quantitative predictions of all of the above-mentioned scattering observables 
at energies lower than 200 MeV. Indeed the latter conclusion is consistent with corresponding 



nonrelativistic microscopic studies 



261 ] . Guided by our phenomenological result, the 



next phase will be to incorporate the PB and density-dependent corrections within the 
context of microscopic relativistic dynamical models in a manner similar to the successful 
nonrelativistic g- folding model developed by the Melbourne group j^ . 

In closing, we emphasize that not only is the RIA a highly predictive model, but it 
also provides a microscopic meson-exchange picture for understanding the behavior of the 
successful global Dirac optical potentials. Another attractive feature of the RIA is the 
consistent application of the relativistic meson-exchange NN scattering matrix for generating 
microscopic folding optical potentials - required for calculating the scattering wave functions 
necessary for evaluating relativistic transition matrix elements - and also for describing the 
driving NN reaction mechanism in relativistic distorted wave models of nucleon-induced 
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reactions. One of the most useful applications envisaged for the RIA will be to generate 
microscopic optical potentials for studying elastic and inelastic scattering of nucleons from 
unstable neutron- and proton-rich nuclei for which global optical potentials do not exist. 
Indeed, existing global Dirac optical potentials have been constrained to reproduce elastic 
proton scattering from stable nuclei 2j, and hence there is no reason to believe that these 
potentials can be reliably extrapolated for studies of exotic nuclei. On the other hand, one 
can readily extend the RIA folding procedure to calculate microscopic optical potentials for 
exotic nuclei. The latter can easily be realized due to the current availability of suitable 
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as well sophisticated relativistic mean field 
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5l|. Future work will focus on studying 



relativistic meson-exchange models j^, Q 
nuclear structure models for unstable nuclei 
elastic proton scattering from exotic nuclei. 
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FIG. 2: (Color online) Differential cross sections dcr/cKl (in units of mb/sr), analyzing powers Ay 
and spin rotation functions Q, plotted as a function of center-of-mass scattering angles 6 c .m. ( m 
degrees), for elastic proton scattering from 208 Pb at incident laboratory kinetic energy of 200 MeV. 



The experimental data 



25[ are denoted by open circles. Uncorrected RIA predictions are repre- 



sented by black dashed curves (legend: RIA; color red online), the black dotted curves correspond 
to predictions based on the EDAI Dirac global optical potentials (GOP) [2| (legend: GOP; color 
olive online), and RIA calculations including phepomenological Pauli blocking (PB) corrections are 
indicated by the black solid curves (legend: PB). 
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FIG. 3: (Color online) Same as Fig. [2l 
98 MeV. The experimental data at 121 



except for incident laboratory kinetic energies of 121 or 



2J and 98 



22| are denoted by open circles. 
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FIG. 4: (Color online) Same as Fig. [21 except for an incident laboratory kinetic energy of 65 MeV. 
The gray dashed curve (color cyan online) corresponds to RIA predictions based on the renormal- 
ized medium- modified <rN and wN coupling constants (denoted by the legend: g*), and the gray 
solid lines (color cyan online) denote the combined effect of medium modified coupling constants 



211 ] are denoted by 



together with PB corrections (legend: g* with PB). The experimental data 
open circles. The analyzing powers A y and spin rotation functions Q are multiplied by 3 at the 
c .m. angles from to 25 degrees. 
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FIG. 5: (Color online) Linear interpolation and extrapolation of PB correction factors calculated 
from a relativistic Dirac-Brueckner approach. The four curves from top to bottom indicate the PB 
correction factors corresponding to the imaginary (Im) and real (Re) parts of the vector V and 
scalar S optical potentials. In the code corresponding to Ref. [8J, the values at Ti a b = 170 MeV 
are chosen same as T\ a \> = 200 MeV, which causes the flat section in the curves. The linear 
extrapolation at T[ a b = 98 MeV is emphasized. 
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FIG. 6: (Color online) (a) Analyzing power A y for elastic proton scattering from Pb at 98 MeV, 
plotted as a function of center-of-mass scattering angles 9 c .m. ( m degrees). The experimental 



data are taken from Ref. 



22j]. Uncorrected RIA predictions are represented by black dashed 



ive online) correspond 



curves (legend: RIA; color red online), the black dotted curves (color ol 
to predictions based on the EDAI Dirac global optical potentials (GOP) [2| (legend: GOP), and 
RIA calculations including extrapolated Dirac-Brueckner and phenomenological PB corrections 
are denoted by gray [legend: PB (DB); color cyan online] and black solid curves [legend: PB (this 
work)], respectively, (b) The corresponding real (Re) and imaginary (Im) scalar S and vector V 
optical potentials (in units of MeV) for elastic proton scattering from 208 Pb at 98 MeV, plotted 
as a function of the nuclear radius r (in units of fm). In particular, the two black solid curves 
corresponding to the real parts of the scalar and vector optical potentials are covered by the gray 



curves. 
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FIG. 7: (Color online) Same as Fig. [2j except for elastic proton scattering from 40 Ca at incident 



laboratory kinetic energy of 152 MeV. The experimental data at 152 MeV are taken from Ref. 



29|]. 
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FIG. 8: (Color online) Strengths of the real (Re) and imaginary (Im) scalar S and vector V optical 
potentials (in units of MeV) at the center (radius r = 0) of a 208 Pb nucleus as a function of incident 
laboratory kinetic energies Ti a b (in units of MeV). The identification of the curves is the same as 
in Fig. HI 
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FIG. 9: (Color online) Same as Fig. 0J except for an incident laboratory kinetic energy of 30 MeV 
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FIG. 10: (Color online) The total reaction cross section 07? (in units of b) for elastic proton 



scattering from 208 Pb versus the incident laboratory kinetic energy T\ a \, (in units of MeV). 
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